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1. Chiral symmetry 

About a decade ago Gade 1 ' investigated the effect 
of chiral symmetry in a nonlinear sigma model (Gade 
model). She showed that this symmetry causes the lo- 
calisation length and the density of state to diverge at the 
band center in 2D. This behaviour at the band center is 
different from that found in the "standard" universality 
classes of Anderson localisation. 

If a Hamiltonian H anti-commutes with a unitary ma- 
trix U then H has chiral symmetry. In this case energy 
eigenvalues occur in pairs of opposite sign with pairs of 
eigenstates ip n and Utp n . When there is no diagonal dis- 
order the Hamiltonian may have chiral symmetry. This 
occurs if we can divide the system into two sublattices 
such that nearest neighbour sites belong to different sub- 
lattices. If this is possible, the transformation U which 
changes the sign of the wave function only on one sub- 
lattice anti-commutes with the Hamiltonian. If the total 
number of lattice sites is odd, chiral symmetry guaran- 
tees the existence of zero energy eigenstate(s), or zero 
mode(s). Such states, if they exist, have special proper- 
ties e.g. they are supported only on one sublattice. The 
precise number of zero modes depends on the difference 
of the number of sites in each sublattice. 2 - ) 

One example of a system with chiral symmetry is the 
random magnetic flux model. Anomalous behaviour is 
observed at the band center. 3 ) The results indicate that 
the entire band is localised except for a critical state at 
the band center. 



»-j ■ 2. The SU(2) model 



Here we report a numerical investigation of localisation 
in the SU(2) model 4 ) without diagonal disorder. The 
Hamiltonian of the SU(2) model has time reversal sym- 
metry but spin rotation symmetry is broken by a random 
spin-orbit coupling. Thus we expect that states will be- 
long either to the standard symplectic class or the chiral 
symplectic class. This model is different from the ran- 
dom magnetic flux model in that the zero mode is sur- 
rounded by a metallic phase, not by an insulating phase. 

The Hamiltonian of the SU(2) model describes non- 
interacting electrons on a simple square lattice with near- 



est neighbour SU(2) random hopping 

(i,j),o;<r' 



(1) 



where ct(cj CT ) denotes the creation (annihilation) oper- 
ator of an electron at site i with spin a. We distribute 
hopping matrices randomly and independently with uni- 
form probability on the group SU(2): 



U(i,j) 



— e -17y sin 8ij e~ lolij cos Oij 



(2) 



where a and 7 are uniformly distributed in the range 
[0, 27t), and 9 distributed in the range [0, n/2] according 
to the probability density, P(9)d8 = sin(20)d0. 

The SU(2) model has a mobility edge at E c ~ 3.253. 4 > 
In what follows we focus attention on the properties of 
the model near the band center at E = 0. We calculate 
the localisation length on a quasi-lD strip and attempt 
to extrapolate to the 2D limit. 

3. Behaviour of the localisation length 

We consider a quasi- ID strip whose width is L and 
calculate the localization length A at arbitrary energy 
E with the transfer matrix method. 5 - ) In the trans- 
verse direction we impose either fixed boundary condi- 
tions (FBC) or periodic boundary conditions (PBC). The 
Hamiltonian has chiral symmetry except when PBC are 
imposed on a system whose width is odd. In this latter 
case chiral symmetry is broken. We analyse the depen- 
dence of the re-normalized localization length A = X/L 
on L. 

The re-normalized localization length A at E = 0.0 
and E = 0.5 as a function of L are shown in Figures 1 
and 2 respectively. The width L ranges from 7 to 128. 
We performed 10 6 ~ 10 8 transfer matrix multiplications 
to obtain data with an accuracy of 0.05% to 2%. 

Figure 1 shows the results at the band center E = 0.0. 
(Note that for odd L and FBC the localisation length 
A diverges 6 - 1 so no data is presented for this case.) A 
striking dependence on the parity of L is observed. As a 
function of L, A approaches a strongly parity dependent 
constant value for L — > 00. For even L this asymptotic 
value is boundary condition independent. To analyse the 
data in detail, we fit the data to 



A = a + cL 



(3) 
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Fig. 1. A vs L at E = 0.0. The solid lines are the fit. 



estimated (asymptotic) slopes for both curves in Figure 2 
are in agreement. 

To compare this behaviour with that in the metallic 
phase of the standard symplectic class, we also performed 
simulations of the SU(2) model with diagonal disorder by 
adding a term 

H d = zA* c ™ (6) 

i.fT 

to the Hamiltonian (1). Here distributed randomly 



Table I. The best fit estimates of fitting parameters at E = 0.0 
with their 95% confidence intervals. Q is the goodness of fit 
probability. 



The results are tabulated in Table I. For even L the 
estimates of the asymptotic value a are the same within 
numerical accuracy. For odd L and PBC the asymptotic 
value is much larger than for even L. The L independent 
behaviour of A at large L is typical of a critical point and 
may indicate that band center is critical in this system. 
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Fig. 2. A vs L at E = 0.5. The solid lines are the fit. 

In contrast at E = 0.5 (Figure 2) the parity de- 
pendence is negligible (except at small L under PBC). 
Though a much strong boundary condition dependence 
is now observed, it seems this will disappear in the limit 
of large L. For PBC the data can be fitted by 

K = a + b\nL. (4) 

and for FBC by 

k = a + b\nL + cL- y . (5) 

The results are tabulated in Table II. In the limit of 
sufficiently large L it seems that a boundary condition 
independent logarithmic increase of A will be recovered. 
Note in particular, by reference to Table II, that the 



Q a c y 

FBC, Even 08 1.39 ± .01 2.6 ± .4 1.12 ± .08 

PBC,Even 0.8 1.381 ± .002 -2 (-6, -0.8) 2.2 ± .5 

PBC, Odd 0.8 4.48 ±.03 3.2 ±1.3 1.2 ± .2 



Table II. The best fit estimates of fitting parameters at E = 0.5 
with their 95% confidence intervals. Q is the goodness of fit 
probability. In fitting we eliminated the data for L < 8. 





Q 


a 


b 


c y 


PBC 


0.1 


4.14 ± .04 


0.653 ± .013 




FBC 


0.2 


5.0 ± 1.6 


0.6 ±.3 


12 ±6 1.0 ±.3 



and independently with uniform probability on the in- 
terval [—1,1]. The diagonal disorder breaks chiral sym- 
metry and the system belongs to the standard symplectic 
class. The re-normalized localization length A is again 
a logarithmic function of L and we are able to fit the 
result by (4) . We also get an estimate of the coefficient 
b = 0.64 ± .02 in the SU(2) model which is consistent 
with those in Table II . This indicates that the states 
away from the band center in the SU(2) model without 
diagonal disorder are metallic and belong to the standard 
symplectic universality class. 
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